arXiv:1508.00720v3 [cond-mat.stat-mech] 30 Sep 2015 


Exact work statistics of quantum quenches in the anisotropic XY model 
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We derive exact analytic expressions for the average work done and work fluctuations in instanta¬ 
neous quenches of the ground and thermal states of a one-dimensional anisotropic XY model. The 
average work and a quantum fluctuation relation is used to determine the amount of irreversible 
entropy produced during the quench, eventually revealing how the closing of the excitation gap leads 
to increased dissipated work. The work fluctuation is calculated and shown to exhibit non-analytic 
behavior as the pre-quench anisotropy parameter and transverse field are tuned across quantum 
critical points. Exact compact formulas for the average work and work fluctuation in ground state 
quenches of the transverse field Ising model allow us to calculate the first singular field derivative 
at the critical field values. 
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I. INTRODUCTION 

A quantum quench (or non-adiabatic change in the 
system Hamiltonian) of an isolated many-body system 
generally leads to a state that is far from equilibrium. In 
particular, a thermal state po = / tr does not 

generally evolve into a canonically distributed (Gibbs) 
state by unitary dynamics under a different Hamiltonian 
Hi m. There is much interest in the study of these highly 
excited quenched states, especially when the quench is 
done about a quantum critical point, because of the pos¬ 
sibility of observing universal phenomena HH, dynami¬ 
cal phase transitions [7] , quantum revivals or thermaliza- 
tion ISHU], and singular behavior in the thermodynamic 
limit [T51 - BD] . Also, when the quench involves interac¬ 
tion parameters that couple different subsystems within 
a system, non-trivial correlations and entanglement may 
be generated between these subsystems pTH23] . Thus, it 
has been suggested that quench-based protocols may be 
used to measure entanglement entropies 

Furthermore, when a quantum quench is viewed as 
a thermodynamic process, it can lead to dynamically- 
generated quantum fluctuations that are described by 
generalized fluctuation relations [251127j . For example, 
the work done on a quenched isolated system is described 
by a probability distribution whose characteristic func¬ 
tion takes the form of a two-time correlation function 
[28j . This correlation function can be expressed in terms 
of free energy differences of a system in effective thermal 
equilibrium uniisn]. This emerging thermodynamic de¬ 
scription therefore opens up the possibility of the use of 
thermodynamic (statistical) methods in the description 
of isolated quantum systems driven away from equilib¬ 
rium (cf. [M [T51 [T71 [T51 [5TH5I] ). 

In this paper, we continue previous investigations on 
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the work done in quenched quantum chains. In par¬ 
ticular, our emphasis is on the theoretical analysis of 
the statistics of work done in arbitrary quenches of the 
ground and thermal states of the XY model following 
the formalism developed in Ref. [28]. Work statistics 
have been studied in other quenched systems such as the 
transverse field Ising model [71|T5J[THl|35j, the XXZ model 
[551157] and integrable field theories |5J [2^ , Luttinger liq¬ 
uids [35] , low-dimensional quantum gases [U [TU |36l EH] , 
critical spin-boson models |31j . and optomechanical sys¬ 
tems [32]. These theoretical studies have recently been 
complemented by experimental interferometric measure¬ 
ments of the probability distribution of work done on 
quenched gases of two-level systems [3DHi2] . Such mea¬ 
surements have also been viewed from a quantum infor¬ 
mation perspective, where the measurement of work done 
is considered as a generalized quantum measurement (a 
positive operator valued measure or POVM) [33j . 


The results of this study can be divided into two main 
parts. In the first part, we report a calculation of the 
average work done in an isolated quenched XY chain 
(Section HI|. When the system is initially in a canon¬ 
ically distributed state po, the quench is associated with 
the production of a so-called irreversible entropy A^irr, 
which may be interpreted as a measure of the irreversibil¬ 
ity of the quench [T31I33H35] . Indeed, in an instantaneous 
quench AS'irr can be expressed mathematically as the 
relative entropy (Kullback-Leibler divergence) between 
the probability distributions of work done in the for¬ 
ward quench protocol and the reversed protocol (with 
initial and final Hamiltonians swapped) [I8j . Addition¬ 
ally, ASiir has been recently proposed as an experimen¬ 
tally accessible quantity that may be used to establish 
the arrow of time in isolated quantum systems [32]. We 
give an exact solution for AS'iri. for arbitrary sets of pre- 
and post-quench Hamiltonian parameters and examine 
sharp increases in irreversible entropy production about 
quantum criticality. 


Second, we calculate the fluctuation in the work done 
in an arbitrary ground state quench and obtain exact con¬ 
tour integral representations in the thermodynamic limit 
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FIG. 1. (Color online) Phase diagram of the XY model show¬ 
ing the quantum critical lines (red solid lines). The anisotropy 
parameter is 7 and the transverse magnetic field is h. In the 
FX and FY phases the ground state has ferromagnetic order 
in the x- and y-directions, respectively. The paramagnetic 
P± ground state has no long-range order on the longitudinal 
a;i/-plane. 


(Section [rv|). We find that the work fluctuation is gener¬ 
ally not an analytic function of the pre-quench (but not 
post-quench) Hamiltonian parameters along the quan¬ 
tum critical lines. For quenches done on thermally mixed 
states, a numerical analysis reveals that this non-analytic 
behavior is weakened in the sense that singular behavior 
begins to appear in higher-order derivatives. These re¬ 
sults suggest that quench protocols may be used to locate 
quantum critical points in systems with unknown phase 
diagrams. 

Past investigations of quenches in the one-dimensional 
XY model have focused on the time evolution (quench 
dynamics) of correlation functions [HOUM] , the Loschmidt 
echo [3], defect production [sniissM], and entanglement 
measures m [HI lini EH]- Detailed theoretical analyses 
of the work statistics have been performed for quenches 
along the transverse field Ising line of the XY model p(5l 
I18| . and we therefore emphasize in this paper the effects 
of anisotropy on the average work and work fluctuations. 
Still, we are able to report some new exact formulas for 
the average work © and work fluctuation ( [M| ) along 
the Ising line to complement these previous studies. The 
overall contribution of our analytic results is therefore 
the completion of the study of work statistics of quantum 
quenches in the full parameter space of the XY model. 


II. MODEL AND DERIVATIONS 


XY model and quench protocol —The one-dimensional 
ferromagnetic XY model Hamiltonian is 


The XY Hamiltonian can be diagonalized by a se¬ 
quence of transforms: a Jordan-Wigner transform to 
spinless fermionic operators, a Fourier transformation 
into momentum space, and a Bogolyubov rotation to ob¬ 
tain a non-interacting fermion theory. The final result 
is 


N-l 


n =0 


( 2 ) 


where xh s-nd Xn are canonical fermion creation and an¬ 
nihilation operators, respectively m- The dispersion re¬ 
lation e{qn) is 

e{qn) = [(/i-cos( 7 „)^-k 7 ^sin^ 9 „]^^^, (3) 

and qn = 2Trn/N for odd N 01 qn = 27r(n -|- ^)/N for 
even N. In the thermodynamic limit V —>■ 00 the effects 
of N being odd or even are negligible. 

The phase diagram of the XY model is given in Fig¬ 
ure The excitation gap vanishes along the critical field 
lines {h = ± 1 ) and the critical isotropic line (7 = 0 and 
\h\ < 1 ) and the system is quantum critical in these re¬ 
gions. 

We consider here quench protocols in which the initial 
Hamiltonian Hq has anisotropy parameter 70 and field 
ho and is instantaneously changed to the final Hamil¬ 
tonian Hi with corresponding parameters 71 and hi. 
This quench is performed with the system initially pre¬ 
pared in a canonically distributed mixed state po = 
e-^^° / ti e-^Ho gx an effective reciprocal temperature 
(3 = {kT)~^, or in the ground state po = |0)(0| of Hq 
at effective zero temperature /3 —> 00 . The system is iso¬ 
lated and is not in contact with a thermal bath during 
the quench. The work done W during the quench is de¬ 
fined as the difference between projective measurements 
of the system energy after and before the quench. 

Work statistics —The statistical properties of the prob¬ 
ability distribution of work done p{W) in an instanta¬ 
neous quench are most accessible from the characteristic 
function 

Giu,P)= j e™’^p(IF) dW = tr po , (4) 

or for a ground state quench G(u) = |0) 

|28) . Using the known eigenstates of the XY model, the 
zero temperature characteristic function G{u) can be cal¬ 
culated and expressed compactly as 

G{u) = ^[cos^ A„ + sin^ A„]. (5) 


N r 


"-si: 

1=1 


1 + 7 


^3 


1-7 


y V 

^3^3+1 


ha^, 


( 1 ) 


with anisotropy parameter 7 , transverse field h, and 
the Pauli a-matrix on site j. We impose periodicity on 
the spin chain so that a%^i = cr“. 


Here, AE = Ei — Eq is the difference between the post- 
and pre-quench ground state energies and A„ = 0i{qn) — 
0oiqn) is the difference between the post- and pre-quench 
Bogolyubov angles that satisfy 


tan[26»i(q„)] 


li sin qn 
hi cos Qn 


( 6 ) 
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The derivation of this characteristic function follows 
closely that in the transverse field Ising case where 70 = 

7i = 1 [T5] . 

The cumulants of the work distribution can therefore 
be generated by repeatedly differentiating the generating 
function In G. In particular, the average work done {W) 
and variance in work done = (IT^) — (IT)^ are 


using the formula kE{k~^) = E{k) — (1 — k?)K{k) |61) . 
Asymptotic series expansions of this result about the pre¬ 
quench critical fields |/io| = 1 reveal that the average 
work per spin is continuous, but not infinitely differen¬ 
tiable, at |/io| = 1- Specifically, logarithmic singularities 
are generally observed in the field derivative about the 
critical field ho = 1: 


iW) = - lim ^InG, 
i u-j-o du 


= — lim 




0 dv ? 


InG. 


(7) 

( 8 ) 


III. AVERAGE WORK DONE AND 
IRREVERSIBLE ENTROPY 

From the characteristic function we can obtain the 
average work done in a finite temperature quench: 


N-l 


= b X! [^o(9n) - ei(9n) cos2A„] tanh 


n =0 


^eo(gn) 


(9) 

In the thermodynamic limit the sum for the average work 
per spin {w) = limAr_>oo(kF)/A^ can be transformed into 
a definite integral that is solvable by quadratures: 






{ho — hi){ho — cos k) 


+ 


£o{k) 

70(70-71) 

eo(fc) 


tanh dk. 


( 10 ) 


For a ground state quench we take /3 ^ 00 before sum¬ 
ming or integrating over wave vectors. This exact expres¬ 
sion reveals that the average work done can be neatly di¬ 
vided into a contribution from the field change {hi — ho) 
and a contribution from the anisotropy change (71 — 70 )- 
Ground state quench along the Ising line —For exam¬ 
ple, in a zero temperature transverse field quench along 
the Ising line 70 =71 = 1, the preceding integral at 
/3 —>■ 00 can be evaluated exactly in terms of complete 
elliptic integrals [60] : 


{w) = 


hp-hi jE{ho)-{l-hl)K{ho), |/io| < 1, 
irho y\ho\E{ho~^), |/io| > 1, 


( 11 ) 

where K{k) and E{k) have the integral representations 


■ ir/2 

(1 - k"^ sin^ du, 

( 12 ) 

■jr/2 

(1 — k"^ sin^ m)^/^ du. 

(13) 


The exact expression (111 for {w) can be cast in the more 
compact form 


{w) = 


hp — hi 
nho 


X |/io| E{ho ^) 


(14) 


d{w) ^ hi-1 ^ I ln(l - ho), ho < 1 , 
dhp 2 tt |ln(/io — 1), 1 ^ hp, 

and about the critical field hp = — 1 : 


(15) 


d{w) _ hi + l fln(l - |/io|), -l<ho, 
dhp 27r ^ |ln(|/io| - 1), hp < -1. 


(16) 


An exceptional case occurs when the post-quench Hamil¬ 
tonian is also at critical field. The non-analyticity at 
ho = /ii = ±1 is weakened in the sense that the deriva¬ 
tive d{w)/dhp exists at critical hp: 


d{w) 

dhp 


ho—±l 


zb — , for hi — zbl, 
TT 


(17) 


and now the logarithmic divergence first appears in the 
second derivative d'^{w)/dhp. To our knowledge, this is 
the first exact characterization of the non-analytic be¬ 
havior of the average work done in a quenched trans¬ 
verse field Ising model at quantum criticality. (We men¬ 
tion, however, that logarithmic singularities in the work 
fluctuation have been reported in field quenches that are 
localized to a single spin site m-) 

Irreversible work and entropy —The average work done 
in a quantum quench also provides information on the ex¬ 
cess work (or irreversible/dissipated work Wirr = (kF) — 
AE) done above the free energy change of an analogous 
reversible process that takes Hp to Hi at reciprocal tem¬ 
perature 13. This free energy change AF is given by the 
Jarzynski relation [29] 


AF =-^lnG{i/3), 


(18) 


and may be expressed in terms of effective canonical par¬ 
tition functions Zi = tre“^^‘ through the Tasaki-Crooks 
fluctuation relation jsniisi]: 


N-l 


AF = - - In ^ V In 

0 Zp 0 ^ ■ 


n =0 


cosh[/3eo('?n)/2] ' 


(19) 


The irreversible work done is positive as excitations, or 
defects, are produced during the quantum quench. It 
is then usually associated with an irreversible entropy 
production according to AS'iir = 0Wi„ = 0{{W) — AF). 

We remark that the expressions derived here for the av¬ 
erage work ([^ and free energy change (191 are identical in 
form to the corresponding equations previously obtained 
for the transverse field Ising model miiiHi. The general¬ 
ization to the full XY model simply involves the modifi¬ 
cation of the dispersion relations and Bogolyubov angles 
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(a) Field quench hi — ho = 0.01 at constant anisotropy 7 . 


(b) Anisotropy quench 71 — 70 = 0.01 at constant field h. 


FIG. 2. (Color online) Irreversible entropy produced per spin Asirr = limjv^oo ASiyy/N at low temperature /? = 100 for (a) 
small field quenches Sh = 0.01, and (b) small anisotropy quenches Sj = 0.01. Black arrows on the field (anisotropy) quench 
graph denote increasing anisotropy I 7 I (field |/i|) magnitudes from 0 to 2.0 in steps of 0 . 2 . 


to account for the anisotropy 7 . With the full expression 
however, we are able to numerically investigate the effects 
of anisotropy on the quenched model and explore critical 
phenomena along the critical XX line. In the examples 
that follow, we study how field quenches are modified by 
anisotropy and how the behavior of anisotropy quenches 
depend on whether the transverse field is above or below 
the critical threshold. 

First, we consider small transverse field quenches 6h = 
hi — ho = 0.01 at fixed anisotropy and low tempera¬ 
ture (3 = 100. Figure shows how the irreversible en¬ 
tropy is pronounced for quenches in the vicinity of the 
critical field due to the increased amount of irreversible 
work done in defect production as the energy gap closes 
[55H57] . Accounting for anisotropy in the XY model 
reveals additional features of AS'irr than previously ob¬ 
served along the Ising line. For instance, the production 
of irreversible entropy is greatest about the multicriti- 
cal points at the intersection of the XX and critical field 
lines (seen as the peaks at /iq = ±1 of the thick solid 
line in Fig. [2^ . Also, we observe the different effect of 
anisotropy on small 5h field quenches deep in the fer¬ 
romagnetic \h\ < 1 and paramagnetic phases \h\ > 1 . 
In the former, where ferromagnetic order is established 
along the x- and y-directions, increasing anisotropy from 
I 7 I = 0 decreases the irreversible entropy production, 
while in the latter increasing anisotropy increases Asin-. 

Second, we take a look at the case of anisotropy 
quenches <57 = 71 — 7 o = 0.01 at fixed field and low 
temperature /3 = 100. In Fig. we see that irreversible 
entropy production peaks when the quench is done about 
the isotropic XX line. However, this peak is most promi¬ 
nent when the quench is done within the critical region 
\h\ < 1, again mirroring the increased dissipated work 
at criticality. However, unlike the previous case of field 
quenches, less irreversible entropy is produced when the 
anisotropy quench is done about the multicritical points 
{h = ±1) than when done across the critical XX line 


at low transverse fields \h\ < 1. That is, in a given 
anisotropy quench about the isotropic XX region, an in¬ 
creasing transverse field suppresses the production of ex¬ 
citations. 

Non-zero temperature —In the limit of infinite temper¬ 
ature /3 —>■ 0 the irreversible entropy AS'irr vanishes triv¬ 
ially as lim^_>,o P{W) and lim^_>,o PAF vanish individu¬ 
ally. In this extreme case, the pre- and post-quench states 
are both maximally mixed and have zero relative entropy. 
In contrast, when the pre-quench state is prepared at 
a finite non-zero temperature, the temperature depen¬ 
dence of the average work ( 10 ) and free energy change 
(19) are difficult to study analytically because of their 
hyperbolic dependence on the energy spectrum cq- Still, 
the integrals for the irreversible entropy and irreversible 
work per spin are readily evaluated by quadrature meth- 



* 0=70 


FIG. 3. (Golor online) Irreversible work per spin wiyy = 
limjv->oo Wiyy/N for small magnitude ground state quenches 
Sh = Sy — 0.01 along the ho = 70 line. Three critical lines are 
crossed, which correspond to peaks in the irreversible work. 
These peaks become sharper as the temperature is decreased 
to zero (the inset shows the peak at /iq = 7 o = 0 ). 
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FIG. 4. (Color online) In small magnitude quenches Sh = 
Sj — 0.01 along the ho = 70 line, the irreversible work peaks 
at some non-zero temperature T* away from criticality (top). 
The black arrow indicates increasing ho = 70 from 0.0 to 1.5. 
The excitation gap A* is of the order of kT* (bottom). 


ods and we performed empirical studies of the tempera¬ 
ture dependence of these quantities. For instance, in the 
case of small magnitude quenches 6h = 67 = 0.01 begin¬ 
ning at ho = 7 o, peaks in the irreversible work per spin 
Wirr = ^Sirr/P Can be examined as three critical lines 
are crossed (F ig. [^ . As in the case of quenches along 
the Ising line |18j these peaks are caused by increased 
defect production as the excitation gap closes at critical¬ 
ity. Additionally, it turns out that the irreversible work 
per spin Wirr saturates to a constant value in the thermo¬ 
dynamic limit at zero temperature, that is, it does not 
diverge at any temperature. This result implies that the 
total irreversible entropy, like a thermodynamic entropy, 
scales extensively at all hnite temperatures ASirr ~ 
as IV —>■ 00 . Furthermore, for these small magnitude 
quenches, the irreversible work vanishes at high helds 
and large anisotropy at all temperatures as the relative 
parameter changes Sh/ho and Jy/yo become smaller |18j . 

Additionally, an early reviewer of this paper has 
pointed out that in Fig. the irreversible work per spin 
Win- does not have a monotonic dependence on tempera¬ 
ture. Since the irreversible work is the difference of two 
quantities (the average work done and effective free en¬ 
ergy change) with different temperature dependence, it 
is reasonable to anticipate such non-monotonic behavior. 
Indeed, observations from several numerical evaluations 
of the irreversible work reveals that Wi„ peaks at some 
non-zero temperature T* away from criticality. For sev¬ 
eral small magnitude quenches, including the previous 


Phase of Ho 

Poles inside C 

Poles on C 

FX, \ho\ < 1, 70 > 0 

0 , 21, 22 

none 

FY, \ho\ < 1, 70 < 0 

0, 23, 24 

none 

P+, ho > 1 

0, 22, 23 

none 

P-, ho < -1 

0, 21, 24 

none 

Grit. 70 = 0 line, |/io| < 1 

0 

22 = 23, 2l = 24 

Grit, field ho = Ij 7 o > 0 

0, 22 

2l = 23 

Grit, field ho = 1, 7 o < 0 

0, 23 

II 

(N 

Grit, field ho = —1, 70 > 0 

0, 21 

II 

Cl 

U? 

Grit, field ho = —1, 70 < 0 

0, 24 

2l = 23 

Multicritical points 

0 

2:1 = 22 = 23 = 24 


TABLE I. Pole structure of the integrand f{z) for the work 
fluctuation per spin cr^ in a ground state quench (21 1 . The 
pre-quench parameters ho and 70 of Hq are finite and the 
post-quench Hamiltonian Hi is not critical. 


case along the ho = 70 line (Fig. |^, we have found that 
the magnitude of the characteristic thermal energy kT* is 
of the order of the energy gap A*. To explain this obser¬ 
vation in physical terms, we propose that small thermal 
fluctuations can lead to increased defect production by 
supplying the system with some initial energy to over¬ 
come the excitation barrier. However, further increases 
in temperature far above T* has the opposite effect until 
Wirr is completely suppressed at infinite temperature, as 
expected. We also observe in Fig. that there is a dis¬ 
continuity in T* that is caused by two competing peaks 
in the graph of Wirr versus temperature. We found it 
difficult to explain this feature, but our numerical results 
suggest that this bimodality is due to the separate contri¬ 
butions of the field change Sh and anisotropy change 
to the average work done (10). That is, in the examples 
that we have studied there is no sharp change in T* as 
a function of the quench parameters when either of the 
two terms in Eq. (10) are absent or neglected. 


IV. WORK FLUCTUATIONS 


When the initial state is the ground state of Hq the 
work fluctuation is 

1 Jv-i 

= - ^ ei(9„)2 sin2 2A„. (20) 

n=0 


In the thermodynamic limit the fluctuation per spin 
= limAT-j-oo T? /N can be represented as a contour 
integral over a positively oriented unit circle C via the 
transformation —>■ z: 


f'^'^biiho - cosfc) - 70 (hi - 

io 

(A ^00), 


cos k)Y sin^k dk 
Stt ’ 

( 21 ) 
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FIG. 5. (Color online) The work fluctuation per spin cP' for 
ground state anisotropy quenches at different fixed transverse 
field {h — 0 ,|, 1 ,|) is not analytic when the pre-quench 
anisotropy is 70 = 0 and \h\ < 1. This behavior reflects 
quantum criticality in the gapless XX region of the XY phase 
diagram. 


with integrand 


/w = 


1 iz^-iru±iz-z±r 


16 

The zeros z± of the numerator are 

'yohi - 7 i/io ± \/(7o^i - 71 ^ 0 )^ - (70 - 71 )^ 


( 22 ) 


Z-^ — 

7o -71 

while the (generally) non-zero poles zi are 

ho + \/ 7 q -f /ig — 1 

7o + 1 

ho - V 10 +hi-1 
7o + 1 

—hp + \/7o + ^0 ~ 1 
7o - 1 

-hp - V 1 I +hi-1 

7o - 1 


zi = 


Z2 = 


Z 3 = 


Zi = 


(23) 

(24) 

(25) 

(26) 
(27) 


It turns out that the contour integral (211 simplifies 
to the sum of at most three residues of / for non-critical 
Hp: 


= I^es/(z). 


(28) 


hd<i 


about the location of non-analyticities of the work fluc¬ 
tuation on the (/iq, 7 o)-plane: The work fluctuation is 
generally not analytic at the points {hp,jp) where Hp is 
critical. An in-depth analysis of the pole structure of the 
integrand f{z) proves this claim. As summarized in Ta- 
blel a pole moves from inside the contour C to outside it 
and another moves from outside the contour to inside it 
as Hp is changed across a critical point. Since the poles 
Zi are distinct away from the phase boundaries, the func¬ 
tion cr^(ft.Oi 7 o) has different functional forms within each 
phase of Hp. The work fluctuation itself is continuous on 
the critical lines of Hp (the poles that swap locations be¬ 
come equal at the critical points), but it is not infinitely 
differentiable along these critical lines. That quantum 
criticality occurs when the poles Zi fall on the unit circle 
is reminiscent of the Lee-Yang mechanism of phase tran¬ 
sitions in the classical Ising model [63] , and the analogies 
between these two phenomena are currently being inves¬ 
tigated. 

We also make the additional observation that field 
quenches restricted to the XX line 70 = 7 i = 0 have 
no work fluctuation (21). In the isotropic XX regime the 


transverse field merely acts as a chemical potential for the 
fermionic quasiparticle excitations of the model. Thus, 
energy eigenstates of an XX Hamiltonian with field hp 
are also eigenstates of any post-quench XX Hamiltonian 
with a different field hi. The resulting probability dis¬ 
tribution for the work done is therefore delta-peaked at 
W — AE = —^N{hi — hp) with no work fluctuation for 
any magnitude of field quench. 

Anisotropy quenches —With an exact formula for the 
work fluctuation at hand, we can consider specific ex¬ 
amples where the presence of non-analyticities reveal the 
criticality of the pre-quench Hamiltonian Hp. For in¬ 
stance, we take the case of anisotropy quenches 70 —>■ 71 
at different fixed fields hp = hi = h. As shown in Fig. 
the work fluctuation has a kink (finite discontinuity in 
the first derivative with respect to 70 ) at the isotropic 
line 7 o = 0 only for field magnitudes below the critical 
threshold |/io| < 1- An explicit calculation for the case 
of quenches toward the Ising line (71 = 1) gives 


da"^ 

dip 


da'^ 

J 

70—>-0+ dlo 

70-s-O- 1 



\h\ < 1 , 
\h\ > 1 , 


(29) 


and we find that the kink only appears along the critical 
XX line. 

Ground state quench along the Ising line —Another ex¬ 
ample in which our solution yields concrete results is 
given by a ground state quench of the transverse field 
Ising model. The work fluctuation per spin has a dis¬ 
continuous field derivative at pre-quench critical fields 
hp = ±I (except when hi = ± 1 ) [M] : 


(t 2 = 


{hi - hp)‘^/8, 


|/io| < 1 , 


(hi - hp)‘^/{8\hp\), |/io| > 1 - 


(30) 


This exact analytic result is one of the main contributions In the exceptional case hi = ±1 the non-analyticity at 
of this paper. It allows us to make a general statement criticality is more subtle. A series expansion of about 
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(b) Quench to (hi,71) = ( 1 , 0 ). 


-2-10 1 2 

h 

(a) Quench to (hi,71) = ( 0 , 1 ). 


FIG. 6. (Color online) Work fluctuation per spin in the 
thermodynamic limit — limjv->oo T?/N for ground state 
quenches that end at (a) the zero-field Ising model, and (b) 
the multicritical point at hi = 1. The work fluctuation is not 
analytic along the critical field lines |ho| = 1 and the critical 
XX line 70 = 0, |ho| < 1. 




ho = ±1 = hi reveals that it is the third derivative 
d^{a‘^)/dhQ that first becomes discontinuous at ho = ± 1 . 
The exact result (30) includes the previously calculated 


approximate asymptotic formula obtained for hi ^ 1 

m 

Quench with fixed post-quench parameters —The non- 
analytic behavior of the work fluctuation when the pre¬ 
quench Hamiltonian is critical is further illustrated in 
examples whe re t he post-quench parameters (hi, 71 ) are 
fixed. In Fig. 6 a we show our calculations for in the 
case of a ground state quench toward the zero-field trans¬ 
verse Ising model, which clearly shows kinks in the work 
fluctuation along the critical field and critical isotropic 
regions. The exact analytic solution used to obtain this 
plot is given in Appendix [A| 

We now discuss in some detail a more involved exam¬ 
ple in which the post-quench Hamiltonian Hi lies on the 
multicritical point (hi, 71 ) = ( 1 , 0 ), so that both Hq and 
Hi are quantum critical. The work fluctuation in this 
case is graphed in Fig. | 6 b[ which reveals kinks along the 
pre-quench critical field ho = —I. The difference in the 
field derivatives of the work fluctuation per spin along 
this critical line is 


da^ 

dho 


ho^-l + 


da"^ 

dho 


/lo-i—1- 


The first derivative of the work fluctuation with respect 
to anisotropy is also discontinuous across the critical XX 
line: 


djo 


7o->-0+ 


djo 


(l-ho)V2, |ho|<l. 


7o->-0- 


0 , 


|ho| > 1, 


(32) 

However, the non-analytic behavior of at ho = 1 is 
more subtle. Inspection of Fig. | 6 b| does not immediately 
reveal singular derivatives at ho = I. In Appendix [B] we 
prove that the first discontinuous derivative of with 


FIG. 7. (Golor online) For quenches along the Ising line 
7o = 7i = 1 toward zero-field hi — 0, an initial non-zero 
temperature smoothens the kink of the work fluctuation per 
spin at the critical pre-quench field ho — 1 (top). However, 
the numerically calculated second derivative of cr^ at non¬ 
zero temperature suggests the existence of a cusp singularity 
at ho = ±1 (bottom). 


respect to the field at ho = 1 is, in fact, the fifth. That 
non-analyticities in cr^ are weakened when ho = hi fall 
on the same critical region can be traced back to the 
contour integral formula (21) and Table [l] We find that 
the pole that contributes to the non-analytic behavior 
at the critical field ho is cancelled by one of the terms 
(z — z±)^ in the numerator of the integrand f{z) when 
hi = ho. 

Non-zero temperature —Extending these results to ini¬ 
tially mixed states at finite fi yields: 



f {hp-hif + (70 - 7 i)^ sin^ k 
[ 2 cosh[/3eo(fc)/2] 

[71 ho - 7ohi -I- (70 - 71) cos sin^ k 


m 


X tanh[/3eo(fc)/2]^ 1^. 


(33) 


We see that singular behavior that can potentially arise 
from the criticality of Hp must come from the second 
term enclosed in large parentheses, which has Cp in the de¬ 
nominator. Near criticality, the temperature-dependent 
factor tanh(/3eo/2) only partially cancels the divergence 
caused by vanishing Cp. We have performed some nu¬ 
merical calculations of this term and found that cr^ is 
generally smooth (with continuous first derivatives) at 
critical Hp. However, as in the cases of weakened non- 
analyticities above, we have found numerical evidence 


































that higher derivatives of cr^ are singular at the critical 
points even at non-zero temperatures. In Figure [7] we 
present the case of a non-zero temperature quench along 
the Ising line with post-quench field hi = 0. The fluctu¬ 
ation per spin is seen to be continuous at = 1, while 
its second derivative with respect to ho diverges. This 
result is noteworthy as the singular behavior of observ¬ 
ables at quantum critical points tend to be washed out 
by thermal fluctuations. 
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V. CONCLUDING REMARKS 


Appendix A: Quench to zero-fleld Ising model 


In this paper, we have studied the statistics of the 
work done in quantum quenches of the one-dimensional 
XY model. We obtained the characteristic function of 
the probability distribution of work for a generic quench 
from initial values ho and 70 to final values hi and 71 
of the magnetic field and anisotropy, respectively. The 
average work done during the quench was calculated 
from this characteristic function and sharp peaks in irre¬ 
versible entropy production were observed in the vicinity 
of the quantum critical lines. In small magnitude field 
quenches, irreversible entropy production was shown to 
be enhanced (reduced) by increased anisotropy in the 
paramagnetic (ferromagnetic) phase. We also derived an 
exact contour integral formula for the work fluctuation 
per spin in ground state quenches and showed how non- 
analyticities can arise at points where the pre-quench XY 
Hamiltonian is critical. Finally, we obtained exact and 
compact solutions for the work statistics in the case of 
ground state quenches along the the transverse field Ising 
line. This latter result allowed us to identify logarithmic 
singularities and finite discontinuities in the field deriva¬ 
tives of the average work and work fluctuation, respec¬ 
tively, at the quantum critical pre-quench fields. 

Note added —Shortly after the submission of this 
manuscript a related paper was published reporting simi¬ 
lar numerical results for finite-length quenched XY chains 
in the presence of three-site interactions |65j . 


When the post-quench parameters are (hi, 71 ) = (0,1), 
which correspond to the zero-field transverse field Ising 
model, the residue formula ([^ gives 


(l-7o)^ 7oho 

FX, 

8(1 + 70 )^ (1+7o)^’ 

1 /8, 

FY, 

1 + 7o(l + 4hQ) + 7o + 7o 

(Al) 

8(1 + 7o)^ 

^ 27 oho[ 2 (I - hl)+-fo - 7 ^] 

P±. 


8(1 + 7o)^a/7o + *0 - 1 


In the paramagnetic region has a Taylor series expan¬ 
sion about 7 o = — 1 and is therefore infinitely differen¬ 
tiable with respect to 70 , that is, 70 = —I is a removable 
singular point. 


Appendix B: Quench to multicritical point 


Let us consider ground state quenches toward the mul¬ 
ticritical point (hi, 71 ) = (1,0). We use the residue for¬ 


mula (28) and expand the work fluctuation about the 
critical field ho = 1 in the ferromagnetic and P-|- param¬ 
agnetic phases of Ho- 


= ^0 


I + 3|7o|-4(ho-l) + 


2(ho - I)^ 

I70I 


+ h»l)U 1 + 3 l>l - 4(A„ -1) + TtyTT _ +,.., 


l7ol 


2|7or 


FX and FY, 


P-f. 


(Bl) 


We conclude that the work fluctuation and its first four 
derivatives with respect to ho are continuous at the crit¬ 
ical field ho = 1. Furthermore, at this critical field the 
fifth derivative has the discontinuity 


d^(a^) 

dh^ 


dHa^) 

dh^o 


15 


ho = i+ 


ho = l~ 


(B2) 
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